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ABSTRACT 



The recent Chandra X-ray observations of the nucleus of M31, combined with earlier 
VLA radio and HST UV spectral measurements, provide the strictest constraints thus 
far available on the nature of accretion onto the supermassive black hole (called M31* 
hereafter) in this region. One of the two newly-detected sources within roughly an 
arcsec of M31* may be its X-ray counterpart. If not, the X-ray flux from the nucleus 
must be even lower than inferred previously. Some uncertainty remains regarding the 
origin of the UV excess from the compact component known as P2. In our earlier 
analysis, we developed a unified picture for the broadband spectrum of this source 
and concluded that M31* could be understood on the basis of the accretion model 
for Sgr A* at the Galactic center, though with several crucial differences. Contrary 
to the 'standard' picture in which the infalling plasma attains temperatures in excess 
of lO^'^ K near the event horizon, the best fit model for M31*, under the assumption 
that the UV radiation is in fact produced by this source, appears to correspond to 
a cool branch solution, arising from strong line cooling inside the capture radius. 
Starting its infall with a temperature of about 10^ K in the post-shock region, the 
plasma cools down efficiently to about 10"^ K toward smaller radii. An important 
prediction of this model is the appearance of a prominent UV spike from hydrogen line 
emission, which for simplicity was handled only crudely in the earlier work. It is our 
purpose here to model this line emission with significantly greater accuracy, using the 
algorithm CLOUDY, and to correctly take into account the attenuation along the line 
of sight. We show that this level of sophistication reproduces the currently available 
multi-wavelength observations very well. Very importantly, we predict a spectrum with 
several additional prominent emission lines that can be used to verify the model with 
future high-resolution observations. A non-detection of the predicted line emission 
from M31* would then tilt the favored accretion picture in the direction of a hot Sgr 
A*-type of model, though with only a single point remaining in the spectrum of M31*, 
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additional observations at other wavelengths would be required to seriously constrain 
this system. 

Subject headings: accretion — black hole physics — Galaxy: center — galaxies: individual 
(M31) — galaxies: nuclei — X-rays: galaxies 
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1. INTRODUCTION 

The nucleus of M31 contains two components, PI & P2, separated by about Q!!5 1.9 pc 
(Lauer et al. 1993). Its optical image is dominated by PI while there is a UV upturn in P2 (King 
et al. 1995). The kinematics of stars near the nucleus (Kormendy & Bender 1999) suggests that 
PI & P2 may be identified as the turning points of an elliptical distribution of stars orbiting the 
central black hole of mass 3.0 x 10^ Mq in a slightly eccentric orbit (Tremaine 1995). However, 
the relatively intense UV emission from P2 remains a puzzle. 

Based on earlier radio observations of M31* (Crane, Dickel, & Cowan 1992), Mclia (1992) 
suggested that its nature may be similar to that of the supermassive black hole, Sgr A*, at the 
Galactic center. Identifying the X-ray object TF56 (Trinchieri &; Fabbiano 1991) with M31*, 
King et al. (1995) proposed a compact nonthermal source at this location. However, subsequent 
high resolution HST observations have partially resolved P2 (Lauer et al. 1998; Brown et al. 
1998), showing that it has a half-power radius of (/.'OG. Lauer et al. (1998) argued that these new 
observations point to the UV upturn as being produced by a tightly bound cluster of early type 
stars. Even so, this picture does not appear to be complete. 

The apparent identification of an X-ray source with M31* (Garcia ct al. 2000) suggested that 
the accretion process in M31* and Sgr A* must be quite different (Liu & Melia 2001). If neither 
of the two sources detected within about an arcsecond of M31* are its counterpart, this would 
place a strict upper limit on its X-ray flux and hence the temperature of the emitting gas. In this 
picture, the UV radiation is an extension of M31*'s radio component, which, however, must turn 
over below ~ 1 keV. 

Instead of the model for Sgr A* in which the accreting plasma reaches temperatures of 10^° 
K or higher (e.g., Shapiro 1973a; MeHa 1992, 1994), M31* appears to correspond well to a second 
branch of solutions in which cooling dominates heating and the gas settles down to a much lower 
temperature during its infall. An optically thin plasma with a temperature of about 10^ — 10^ K, 
a typical value at the capture radius rc = 2GM/v^ of a supermassive black hole accreting from 
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stellar winds (Coker & Melia 2000), lies on the unstable portion of the cooling curve (Gehrels & 
Williams 1993). Here, v^o is the velocity of the ambient gas flowing past the central black hole. 
Liu & Melia (2001) showed that in this case two branches of solutions exist, distinguished by 
the relative importance of cooling versus compressional heating at rc- Depending on the initial 
temperature T[rc) and the mass accretion rate M, the plasma either settles onto a hot branch 
or a cool branch, characterized by temperatures of ~ lO^*' K and ~ 10^ K, respectively, toward 
smaller radii. 

Comparing these two solutions to the observations, Liu &: Melia (2001) argued that M31* 
must lie on the cool branch, and that we should therefore see a spectral signature of its cooling 
flow, particularly a prominent UV spike due to hydrogen line emission and soft X-ray line 
recombination. The purpose of this paper is to present a significantly more detailed calculation 
of the broadband emission spectrum, correcting also for the attenuation along the line of sight. 
In section ^, we discuss the standard Bondi-Hoyle model for supermassive black holes accreting 
from the ambient medium. In section |3| & ^ we discuss the specific details of our line emission 
calculations and we identify several emission lines that may be used to verify our model. Section 
P presents the conclusions. 

2. LOW ANGULAR MOMENTUM ACCRETION FLOW 

In this section, we briefly summarize the Bondi-Hoyle accretion model for a low luminosity 
accreting supermassive black hole embedded within a gaseous environment. The latter may be 
produced by the winds of nearby stars. Gas flowing past the black hole will be captured at a 
rate determined by its velocity and specific angular momentum relative to the compact object. 
Because winds far from the stars are always highly supersonic, many shocks can form where 
the winds collide before the gas is captured within the central gravitational field. These shocks 
dissipate most of the kinetic energy and heat the plasma to a temperature of ~ (3 Voo/A)'^ fi/3Rg, 
where Rg is the gas constant and fi is the molecular weight per particle, which is about 1/2 for a 
fully ionized plasma. During the infall, the plasma accelerates under the influence of gravity and 
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again goes transonic. We base our analysis on the global structure inferred for these flows from 
existing 3D hydrodynamical simulations of the infalling gas (e.g., Coker & Melia 1997), and we 
take our starting point to be a radius where the plasma has already re-crossed into the supersonic 
region. Our analysis will adopt a purely radial flow, for which consistency requires a negligibly 
small specific angular momentum in the accreting gas. 

Due to their low luminosity, supermassive black holes accreting in this fashion can be treated 
with significant simplification. For example, the energy loss due to radiation can be treated with 
just a cooling term in the energy equation and the related effect on the dynamics of the flow is 
negligible. In this sense, the structure of the accretion flow and the radiation may be handled 
separately. 



2.1. Dynamical Equations 

The stress-energy tensor for a radial accretion flow with a frozen-in magnetic field is given by 
(Novikov & Thorne, 1973): 

Stt 

where e = ank},T/2ii + niripCp is the energy density. The general expression for a is given by 
Chandrasekhar (1939): 

"-"^l AK,[x) V+H 4i^2(y) ) ' 

where x = rrieC? /k},T, y = rUpC^ /kbT and Ki is the ith order modified Bessel function. Also, n is 
the baryon number density, and p = nkbT/n is the gas pressure. The other symbols have their 
usual meaning. The shear and bulk viscosity has been neglected in the equation because the flow 
is spherically symmetric. In the case that the radiation is not very strong, the radiation can be 
treated as an external field, which just contributes a cooling term to the accretion flow. So we 
don't include it in the stress-energy tensor either. 
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In Schwarszchild coordinates, the metric is given as 

^ ^^2 = ^) ^2 ^^2 ^ + r\de^ + sin^ d,^^) , (3) 



where M is the mass of the supermassive black hole. Then we have u" = (li*, u^, 0, 0) = 
(dt/dr, dr/cdr, 0,0), which is dimensionless in om notation. Because the magnetic field in 
a spherical accretion flow is dominated by its radial component (Shapiro 1973b), one has 

= (B* , B'' ,0,0). The electric field is always zero in the comoving frame of a fully ionized 
plasma with frozen-in magnetic field, so B is orthogonal to u: 

B"ua = . (4) 

The dynamical equations for the radial accretion flow are as follows: 

(nu");„ = , (5) 
T'^t^.f^ = -F" , (6) 

where F° = T^J^,^ is the divergence of the energy-momentum tensor Tj-ad for the photons, which 
gives the plasma's energy and momentum loss rate due to radiation. The cooling rate in the 
comoving frame of the flow is given by A = —cF'^Ua- 

For a steady flow, the continuity Equation (|5|) gives 

M = —Attt'^c u^n nip , (7) 

where the accretion rate M is a constant. 

The law of local energy conservation is given by projecting Equation onto ti", for which 

cu%".^^ = A , (8) 

that is, 

d /e\ d 1 r-A 

CLT \n J drn n 



where 



r = C ^i^T^ag = "4^(2 B^Ur + r Ur B B') (10) 
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is the heating term due to the annihilation of the magnetic field. As usual, we will use prime to 
denote a derivative with respective to r and an over-dot to denote a derivative with respective 
to t. T^fg = {2B^u°'u'^ + B'^g^'f^ - 2S"5^)/87r is the Maxwell stress-energy associated with the 
frozen-in magnetic field. 

The Euler equations are given by 

{g""^ + u''u^)Tf3'^.^ = -F"" + ^u'' . (11) 
Only the radial component of the equation is non-trivial. It gives 

(p + e)u:^u" = -{g^P + u^u^)p,f, , (12) 

where we have neglected the radiation pressure and magnetic pressure, which are much smaller 
than the thermal pressure. 

Then the equations for the radial velocity and temperature are: 

^ d ^ _ _ p' f , 2 GM \ GM 

T' r-A I (u' 2\ 
a— = — — — + - • (14) 

T Kb T n C Ur fl \Ur T ) 

Because the spacetime geometry is dominated by the supermassive black hole, all of these 
equations ignore the effects of self-gravity within the accreting gas. These equations are readily 
solved numerically to provide the radial structure of the infalling gas once the boundary conditions 
and the magnetic field are specified. 



2.2. Boundary Conditions and Magnetic Field 

In principle, a precise knowledge of the wind sources surrounding the accretor would provide 
all the necessary conditions at the capture radius to determine the structure of the infalling gas 
uniquely. However, only some global properties of the ambient medium are known adequately well. 
For example, at the Galactic center, the average density and flow velocity have been determined 
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from emission-line measurements, but even there, a lack of information on the de-projected 
distance along the line of sight to the accretor prevents us from positioning the sources with 
certainty around the black hole. The ambient conditions in the much more distant nucleus of 
M31 are even more difficult to identify with precision. Our approach will instead be to use the 
accurately known spectrum produced by the accreting gas to guide our selection of the boundary 
conditions, which therefore become adjustable parameters for the fit. 

Taking our cue from the conditions prevalent at the Galactic center, we adopt a wind velocity 
Voo of about 500 km s~^, which then also determines the downstream plasma temperature for a 
strong shock: 

T = (3 i;oo/4) • (15) 

Typically, T ^ 3 x 10^ K. Three-dimensional hydrodynamic simulations show that the region 
between the shocks and the transonic point is more or less isothermal. So we would expect the 
temperature at the outer boundary of our modeled region to be about several million degrees 
Kelvin. 

HST observations (Lauer et al. 1998) show that the extended UV source in P2 has a radius 
of about 0.2 pc (~ lO^r^, where rs = 2GM/c^ is the Schwarszchild radius of the central black 
hole). The outer boundary should have roughly the same scale as this. To ensure that this outer 
boundary is located below the transonic point (since we are limiting our domain of solution to this 
region — see above), the radial velocity of the flow must be larger than its thermal velocity at that 
point. We have already limited the temperature at the outer boundary. This therefore provides 
a lower limit for the radial velocity. At the same time, the flow must be bounded, so its radial 
velocity cannot exceed the free fall value at this radius. 

With M and Ur known, the accretion rate then depends on the ambient gas density. Although 
Ciardullo et al. (1998) have argued that the ionized gas density in the nucleus of M31 changes 
dramatically, dropping from > 10^ electrons cm^'^ at a radius of ~ 7" or 23 pc to 10^ electrons 
cm-3 at 1' or 200 pc, the actual density near the capture radius may still be quite different 
from this given the differences in length scales. As such, the accretion rate, or the gas density at 
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the outer boundary, remain as adjustable parameters. 

The structure and intensity of the magnetic field are also poorly known in this context, 
given the uncertain nature of magnetic field dissipation within a converging fiow. Following the 
assessments of Kowalenko & Melia (1999), and their application to Sgr A* in the outer region of 
the accretion flow (Coker & Melia 1997), we will invoke a sub-equipartition field compared to the 
thermal energy density. That is, we will take S^/Svr = fS^a n ki, T/2fi, where < 1. 

2.3. Dichotomy of the Accretion Profiles 

Figure [l| shows the cooling curve for an optically thin plasma with cosmic abundances. The 
thin solid curve is that given by Gehrels & Williams (1993), which we used in our previous 
calculation. The thick solid curve gives a more accurate cooling curve based on more detailed 
modeling with CLOUDY (Ferland 1996). This is the curve used for the calculations reported in 
this paper. It is evident that the plasma is located in an unstable portion of this curve when its 
temperature lies between 10^ and 10^ K. As we have seen above, the temperature downstream of 
the capture radius is typically several million Kelvin. Thus, depending on the gas density in this 
region, which determines the cooling rate at the outer boundary and the accretion rate for a given 
radial velocity, the plasma may either heat up to about 10^*^ K, or cool down to around 10^ K. 

Figure ^ illustrates this dichotomy. The thick solid curve corresponds to our best fit model 
for the multi- wavelength observations of M31*, which has an outer boundary Vq of lO^rg (~ (X'OS), 
a radial velocity of one half of the free fall velocity, an accretion rate of 1.5 x 10^^ g s~^ and a 
magnetic parameter = 0.003. The outer boundary temperature is given by the ratio of the 
thermal to gravitational energy density, which is 0.068 for the best fit model. To emphasize the 
accretion rate dependence of the solutions, we also show in this figure the temperature profiles for 
two configurations with different accretion rates, which are 1.5 x 10^^ g s~^ and 1.5 x 10^'^ g s^^ 
for the dashed curve and the thin solid curve, respectively. Clearly, the solutions depend rather 
sensitively on the accretion rate. We will return to this point in section |5| below. It appears that 
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in order for M31* to simultaneously account for both the radio and UV spectral components, the 
post-shock gas initiating its infall toward the central accretor is highly unstable to cooling. For 
this reason, the temperature of the plasma in the best-fit model drops quickly from ~ 10^ K in the 
post-shock/capture region to ~ 10^ K for the remainder of its inward flow. 

In the hot branch, the cooling is dominated by thermal bremsstrahlung radiation. We would 
expect a flat spectrum which cuts off at Uc ~ kf^ T/h. Figure ^ shows the spectra corresponding 
to the three configurations in Figure ^. In this figure, the two more constraining optical upper 
limits are from Lauer et al. (1998) and the three near Infrared upper limits are from Corbin 
et al (2001). The other data are discussed in Liu & Melia (2001). Note that these spectra do 
not include the contribution from cyclo-synchrotron emission, which is absent in the cool branch 
solutions, but would produce an additional peak at radio frequencies for the hot branch profiles. 
The magnitude and location of this cyclo-synchrotron peak would depend on the strength of the 
magnetic field. Clearly, the Chandra X-ray upper limit (or detection, depending on whether or 
not one of the sources near M31* is its counterpart) already rules out the hot branch solution for 
M31*, under the assumption that the radio component extends into the UV region, even without 
the additional problems introduced at radio wavelengths by this additional component. The cool 
branch solution, on the other hand, accounts very well for the radio and UV emission, and it 
features strong line emission due to the dominant line cooling component in the cooling function. 
Even so, the configuration with the low accretion rate can't produce the observed radio and UV 
flux, so there is a restricted range of parameter values at the outer boundary that appear to be 
consistent with the data. 

3. CALCULATING THE EMISSION SPECTRUM 

The emission spectrum of a cooling gas is a function of the density n and temperature T, 
which in turn are both functions of radius. With the assumption of spherical symmetry, the mean 
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emissivity (in units of energy per unit time per unit volume per unit frequency interval) is 

47ri,(r,^) = A7re{r,e) ij^''^{n{r),T{r)) + ^ /-(n(r), r(r)) 5{,y - iyuUr,9))\ , (16) 

V. lines ) 



where 



^ ' ' [{l + v{r)cose/c) ^ ^ 



is the special relativistic correction to the monochromatic flux of moving emitters (Corbin 



19971 ); 47r jj?'"^™(n,r) ~ 6.8 x 10"^*^ ey^p (-hu/hT) erg s"^ cm'^ Hz'^ is the mean 



emissivity due to thermal bremsstrahlung ( Frank, King, fc Raine 1992 ), with the ionization 



fraction A = (1.0 + 0.445 T~^e^^^'^^''/"^)~^ obtained following the prescription given in Rossi et al. 
(1997); y™'^(n,T) is the mean emissivity (in units of energy per unit time per unit volume) due to 
line emission; and 

g / 1 — v{r) cos ( 



'^iUr,o)=4nA \r ] (18) 



is the Doppler-corrected line frequency, with u^^^^ being the rest-frame line frequency. General 
relativistic effects have been ignored as the vast majority of emission occurs at radii r ^ rs- 

Computationally, we view the inflow as a series of concentric spherical shells, each of which 
is characterized by a set of physical parameters. The velocity v{r) and density n{r) are specified 
by the equations of continuity and motion, and the temperature T(r) is specified by the energy 
balance equation, including terms for line and bremsstrahlung cooling ( Pehrels fc Williams 199*3 ) 



The central mass is assumed to be M = 3 x 10'' Mq and the mass accretion rate is set at 
M = 1.5 X lO^^g s~^, corresponding to the best fit model temperature profile shown in Figure ^. 

We model each shell as a coronal equilibrium zone (appropriate for gas which is mainly 
collisionally ionized) and use the photoionization code CLOUDY ( [Ferland 1996 ; CLOUDY 94.00) 



to compute the line emissivities. Cosmic abundances are assumed throughout, with no grains 
present due to the high temperatures. Radiative transfer effects between the zones are neglected, 
as gas is optically thin above O.lGHz. 

To obtain a spectrum of the flux density, the emissivity is integrated over the entire flow 
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volume: 



= ^ — / 4:7rji, sin 9 drdOdd) , 



(19) 



where Dm31 = 784 kpc is the distance to M31 (Stanek & Garnavich 1998). The flux density is then 
corrected for extinction along the line of sight due to absorption ( [Cruddace et al. 1974 ; Morrison 
fc McCammon 1983D and reddening ( pardelli, Clayton, &: Mathis 1989]) , assuming a column depth 
of A^H = 2.8 X 10 cm" based on Chandra observations (Garcia et al. 200q ), cosmic elemental 
abundances ( palgarno fc Layzer 198^ ), Ay = 0.24 ( |Lauer et al. 1993| ), and E{B - V) = 0.11 
(IKing, Stanford, fc Crane 199E| ). 



4. RESULTS 

Figure ^ shows the intrinsic flux spectrum calculated with our model. The spectrum (and 
cooling) is dominated by the UV recombination lines of warm Fe, O, Si, and S species present in 
the outer flow region. 

Figure |5| shows the model spectrum corrected for extinction. It is evident that the strongest 
emission lines are effectively extinguished by the intervening column of gas. The surviving 
optical-UV and soft X-ray components both show very strong line emission superimposed on a 
bremsstrahlung continuum, and should be easy to distinguish from, e.g., the composite spectrum 
of a stellar population. For comparison, the spectra of T = 10, 000 K and T = 30, 000 K stars 



computed from Atlas stellar atmosphere modeling ( Kurucz 1991 ) are also presented in the flgure. 

An expanded view of the modeled UV and X-ray spectra are shown in Figures ^ and |^. High 
spatial resolution observations of M31* at these wavelengths may provide a very straight-forward 
test of this cooling accretion flow picture. In particular, the UV spectrum is characterized by very 
strong line emission; Table H lists the intensities of the strongest emission lines relative to Lya. 
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5. DISCUSSION 



As mentioned in section these solutions are rather sensitive to the outer boundary 
conditions, particularly the accretion rate. For these simulations, we have used a simplified scheme 
in which compressional and magnetic heating balance line cooling in the outer region. In reality, 
turbulent heating may also play an important role in stabilizing the temperature. Additional 
heating may be provided by the disturbance generated by the stellar population around P2. In 
any event, the essential feature of this model is that regardless of what happens at large radii, line 
cooling can dominate over the heating processes at small radii when the accretion rate is large 
enough. This inner cool region is important because the observed weak radio emission does not 
allow the plasma to be heated up. The two radio measurements of M31* give a flux density of 
28 — 39 /uJy at 3.6cm (Crane et al. 1992 & 1993). The 3.6cm radio flux from an optically thick 
spherical emitter can be estimated using the following equation: 

(3.6 cm) 



-Da/31 rsJ A2 



Thus, a source with radius 10 rg and a temperature of several times lO^'' K will produce far too 
much flux at 3.6 cm. On the other hand, a hot optically thin source cannot produce the strong UV 
emission via self-Comptonization of the radio photons. If the UV component is in fact associated 
with M31*, this argument also rules out many other models developed in recent years to account 
for the emissivity of supermassive black holes in the cores of nearby galaxies, including Sgr A* 
at the Galactic center (the ADAF model, Narayan, Yi, &: Mahadevan 1995; the ADIOS model, 
Blanford & Begelman 1999; and the jet model, Falcke & Markoff 2000). However, if the UV 
emission is not associated with M31*, a hot optically thin plasma may still account for the radio 
emission, as it does for Sgr A*. Distinguishing between these two scenarios via the observation of 
UV lines in M31*'s spectrum makes our proposed measurements highly desirable. 

In calculating the spectrum for the hot branch solution in Figure |^, we have not included the 
effects due to a magnetic dynamo, should the gas circularize before it reaches the event horizon. 
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In the case of Sgr A*, the dominant contribution to the mm and sub-mm spectrum is made by a 
Keplerian structure within the inner 5 rg or so (Meha, Liu & Coker 2000, 2001). In the case of 
M31*, this region may or may not exist, depending on how much specific angular momentum is 
carried inward. For the same conditions as in Sgr A*, M31* would not have a Keplerian region 
since rs(M31*) ~ 10 r5'(SgrA*), i.e., the gas would not circularize before crossing the event 
horizon. So the model for Sgr A* may not work here. The much larger accretion rate implied by 
our model also rules out the possibility of a standard optically thick disk being present in M31*. 
As shown in Figure ^, even with an accretion rate more than 6 orders of magnitude smaller than 
the accretion rate we infer for the best fit case, a small optically thick disk still produces too much 
optical radiation compared to the observations. Such an optically thick disk would not be able 
to account for the UV upturn either, making it unlikely that such a structure is present in the 
nucleus of M31. 

In Figures 3 and 8, we have treated the Chandra flux measurements of the two nearest X-ray 
sources to M31* as upper limits, since neither of these has been established with confidence as 
the actual counterpart to the radio source. A slightly different fit, assuming that the so-called 
"southern" source is indeed the counterpart, was presented in Liu & Melia (2001). None of the 
conclusions presented in this paper are affected qualitatively in this instance, except for slight 
differences in the recombination line intensities. However, if the "northern" source turns out to be 
the actual counterpart (see Garcia et al. 2000), the cold accretion model will not work for M31*. 
A close inspection of Figure 8 shows that the X-ray spectrum (indicated by the butterfly) is then 
too hard to be described adequately by bremsstrahlung emission. The nature of M31* is therefore 
still rather uncertain, making further study of this very interesting object highly desirable. 
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Table 1. Attenuated Emission Line Intensities 



Line^ 


Intensity* 


Lya A1216 


1.00 


VI AA1032, 1038 


0.59 


Mg II AA2796, 2804 


0.47 


CIV AA1548,1551 


0.16 


N V AA1239, 1242 


0.13 


IV A1399 


0.13 


S VI AA933, 945 


0.07 


C 11 A2326 


0.06 


C III A977 


0.06 


Ne VI A999 


0.06 


V A1215 


0.03 


Si III A1207 


0.03 


Si III A1891 


0.03 


Ly/3 A1026 


0.02 


Ly7 A973 


0.02 


C I A1657 


0.02 


N III A991 


0.02 


III AA1661, 1666 


0.02 



'^All wavelengths in A. 
^Relative to Lya A1216. 
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Fig. 1. — Temperature dependence of the cooling coefficient, defined to be the cooUng rate divided 
by the baryon number density squared, for an optically thin plasma with cosmic abundances. The 
thick solid curve is the coefficient used in this paper. The thin solid curve is the result given by 
Gehrels & Williams (1993), which we used in our previous calculation (Liu k, Melia 2001). 
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Fig. 2. — Temperature profile for the accreting gas, for three values of accretion rate. The initial 
temperature in all cases is 1.24x lO^K . The accretion rate are l.Ox lO^^g (thin, solid curve — the 
hot branch solution), 1.0 x 10^'^g s^^ (dotted curve) and 1.0 x lO^^g s~^ (thick, solid curve — the 
best fit model for M31*). Also shown here is the run of proton density (dashed curve) of the best 
fit model, whose scale appears on the right hand side. 
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Fig. 3. — A comparison of the overall spectrum arising from each of the three cases depicted in 
Fig. |2[ The line types have the same definitions as those in the figm'e. The X-ray upper limits 
correspond to the so-called "southern" source detected near M31* by Chandra, whereas the butterfly 
shows the measured spectrum for the "northern" source. 
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Fig. 4. — Modeled intrinsic bremsstrahlung (thin curve) and bremsstrahlung plus line (dark, heavy 
curve) spectra. 
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Fig. 5. — Modeled line plus bremsstrahlung spectrum corrected for extinction and reddening (dark). 
Shown for comparison are stellar atmosphere spectra for T = 10, 000 K and T = 30, 000 K stars 
(thin curves), arbitrarily scaled. See text for description. 
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Fig. 6. — Modeled ultraviolet spectrum. 



-26- 




E (keV) 



Fig. 7. — Modeled X-ray spectrum, rebinned in equal logarithmic intervals. 
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log,o / Hz) 

Fig. 8. — Spectra of optically thick disks. The disks have an inclination angle of 60°, an inner 
boundary of 3rs and the accretion rate is 5.0 x 10^^ g s~^. The thick solid line corresponds to a 
small disk with an outer boundary of 10r5, while the dotted line corresponds to a disk with an outer 
boundary of 1000r5. The thin solid line corresponds to the best fit radial accretion model, though 
here only the hydrogen line emission is included. As in Fig. 3, the X-ray upper limits correspond to 
the so-called "southern" source detected near M31* by Chandra, whereas the butterfly shows the 
measured spectrum for the "northern" source. In this best fit, we have treated these X-ray data as 
upper limits to the actual spectrum of M31* (see text). 



